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FORMATION OF HIGH-REDSHIFT {Z > 6) QUASARS DRIVEN BY NUCLEAR STARBURSTS 
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ABSTRACT 

Based on the physical model of a supermassive black hole (SMBH) growth via gas accretion in a 
circumnuclear disk (CND) proposed by Kawakatu & Wada (2008), we describe the formation of high-z 
[z > 6) quasars (QSOs) whose BH masses are Mbh > IO^Mq. We derive the necessary conditions to 
form QSOs at z > 6 by only gas accretion: (i) A large mass supply with Mgup > IQ^^Mq from host 
galaxies to CNDs, because the final BH mass is only 1 — 10% of the total supplied mass from QSO 
hosts, (ii) High star formation efficiency for a rapid BH growth which is comparable to high-z starburst 
galaxies such as submillimeter galaxies (SMGs). We also find that if the BH growth is limited by the 
Eddington accretion, the final BH mass is greatly suppressed when the period of mass-supply from hosts, 
tsup is shorter than the Eddington timescale. Thus, the super-Eddington growth is required for the QSO 
formation as far as tsup, which is determined by the efficiency of angular momentum transfer, is shorter 
than ^ 10® yr. The evolution of the QSO luminosity depends on the redshift Zj at which accretion onto 
a seed BH is initiated. In other words, the brighter QSOs at z > 6 favor the late growth of SMBHs (i.e., 
Zi K, 10) rather than early growth (i.e., Z{ « 30). For Zi ~ 10, isup — 10® yr is shorter than that of the 
star formation in the CND. Thus, the gas in the CND can accrete onto a BH more efficiently, compared 
with the case for Zi « 30 (or igup ~ 10^ yr). Moreover, we predict the observable properties and the 
evolution of QSOs at z > 6. In a QSO phase, there should exist a stellar rich massive CND, whose gas 
mass is about 10% of the dynamical mass inside ~ 0.1 — Ikpc. On the other hand, in a phase where 
the BH grows (i.e., a proto-QSO phase), the proto-QSO has a gas rich massive CNDs whose gas mass is 
comparable to the dynamical mass. Compared with the observed properties of the distant QSO SDSS 
J1148-I-5251 observed at z — 6.42, we predict that SDSS J1148+5251 corresponds to the scenario of the 
late growth of SMBH with Zi ~ 10, which is accompanied by a massive CNDs with Mg « 5 x IO^^Mq 
and the luminous nuclear starburst LgB at infrared band with Lsb ~ lO^^'ergs"^. Moreover, we predict 
that the progenitor of SDSS J1148-I-5251 can be the super-Eddington object. These predictions can be 
verified by ALMA, SPICA and JWST. 

Subject headings: black hole physics — early universe — galaxiesiactive — galaxies:nuclei — 
ISM:structure — galaxiesistarburst 



1. INTRODUCTION 

Supermassive black holes (SMBHs) with masses in the 
range of 10^ — IO^Mq are the engines that power active 
galactic nuclei (AGNs) and quasars (QSOs). There is also 
ample evidence that SMBHs reside at the center of most 
galaxies (e.g., Kormendy & Richstone 1995; Richstone et 
al. 1998; Ho et al. 1999), including the Milky Way (e.g., 
Genzel et al. 1997; Schodel et al. 2002; Chez et al. 2003). 
In the local Universe, there are tight correlations between 
the masses of SMBHs and the masses and velocity disper- 
sions of the spheroidal components (bulge) of the hosts 
(e.g., Kormendy & Richstone 1995; Richstone et al. 1998; 
Ferrarese & Merritt 2000; Gebhardt et al. 2000; Tremaine 
et al. 2002; Marconi & Hunt 2003; Hiiring & Rix 2004; 
Barth et al. 2005). These observational facts suggest that 
an intriguing link between the formation of bulges and 
SMBHs. 

High redshift QSOs are essential to understand the for- 
mation and evolution of SMBHs. In recent years, more 
than thirty QSOs at z w 6 have been discovered (Fan et 
al. 2000, 2001, 2003, 2004, 2006; Goto 2006; Jiang et 
al. 2006; Willott et al. 2007; Jiang et al. 2009). The 
higher luminosity of the distant QSO at z = 6.42, cor- 



responding to SDSS J1148-h3251 (Fan et al. 2003) with 
~ 10^^ergs~\ implies that a SMBH with mass > lO^M© 
is already in place within ~ 1 Gyr after the big bang, by 
assuming the Eddington luminosity. These requirements 
set significant constraints on the evolution and formation 
of SMBHs in the early Universe. There are some analytical 
and semi-analytical studies (Haiman & Loeb 2001; Haiman 
2004; Yoo & Mirald-Escude 2004; Shapiro 2005; Volonteri 
& Rees 2006; Tanaka & Haiman 2009) and numerical sim- 
ulations (Li et al. 2007; Sijacki et al. 2009) that discuss 
the formation of > IO^Mq SMBH at z > 6. A key physical 
process governing evolution of QSOs is the mass accretion 
toward a SMBH, although BH growth via merging may be 
expected at the high mass BH with A/bh > lO^M© (e.g., 
Shankar 2009 and references therein) . Although the accre- 
tion process is a complicated phenomenon over nine orders 
of magnitude in the size scale, the Eddington-limited ac- 
cretion rate was simply assumed in most previous studies. 
However, the evolution of the SMBH is not only controlled 
by the accretion processes in the vicinity of the SMBH, but 
also related to the mass supply from the host galaxy to the 
circumnuclear region. 

A number of mechanisms for gas accumulation on the 
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~ 1 — 10 kpc scale down to the galactic central region have 
been proposed, e.g., the tidal torque driven by the major 
and minor mergers of galaxies (e.g., Toomre & Toomre 
1972; Mihos & Hernquist 1994, 1996; Saitoh & Wada 2004; 
Saitoh ct al. 2009) and the stellar bars (e.g., Noguchi 
1988; Shlosman et al. 1990; Fukuda et al. 1998; Fukuda 
et al. 2000; Maciejewski et al. 2002; Namekata et al. 
2009), gas drag and dynamical friction in the dense stellar 
cluster (Norman & Scoville 1988) and the radiation drag 
(e.g., Umemura et al. 1997; Umemura 2001; Kawakatu & 
Umemura 2002). However, the accumulated gas does not 
directly accrete onto a SMBH, since the angular momen- 
tum of the gaseous matter cannot be thoroughly removed. 
Thus, some residual angular momentum would terminate 
the radial infall, so the accreted gas forms a reservoir, i.e., a 
circumnuclear disk (CND), in the central ~100 pc around 
a SMBH whose scale depends on the angular momentum 
of the gas. If the gravitational instability takes place in a 
CND, star formation is naturally expected in this region. 
Active nuclear star formation (< 100 pc) has been actually 
observed in nearby AGNs (e.g., Imanishi & Wada 2004; 
Davies et al. 2007; Watabe et al. 2008; Chen et al. 2009; 
Hicks et al. 2009). It is theorized that the nuclear star- 
burst could obscure some types of AGNs (e.g., Ohsuga & 
Umemura 1999; Wada & Norman 2002: hereafter WN02; 
Thompson et al. 2005; Watabe et al. 2005; Ballantyne 
2008; Schartmann et al. 2008; Wada et al. 2009). The 
nuclear starburst also affects the growth of SMBHs, be- 
cause the radiation and/or supernova feedback from star- 
bursts can enhance the mass accretion onto a SMBH (e.g., 
Norman & Scoville 1988; Umemura et al. 1997; WN02; 
Vollmer & Bcckert 2003; VoUmer, Bcckcrt & Davis 2008; 
Collin & Zahn 2008). In order to reveal the final rate of 
mass accretion to the BH region, it is crucial to link mass 
accretion processes from a galactic scale with those from 
an accretion disk in the vicinity of a central BH, via the 
CND. Recently, we have proposed a theoretical model of 
a nuclear starburst disk supported by the turbulent pres- 
sure led by supernova explosions (Kawakatu & Wada 2008: 
hereafter KW08). In KW08, the turbulence excited by su- 
pernovae transports the angular momentum. We showed 
how a SMBH grows from a seed BH, taking into account 
the mutual connection between the mass-supply from a 
host galaxy and the physical states of the CND accompa- 
nied by the star formation. This theoretical model should 
be confirmed for the formation of high-^; {z > 6) QSOs 
whose masses are ~ 10^ Mq. 

This paper is organized as follows. We briefly outline the 
model in KW08 in §2. By adopting KW08, we will show 
the necessary conditions to form QSOs at x > 6 in terms 
of the total accreted gas mass from host galaxies and star 
formation efficiency in the CNDs (§3). In addition, we will 
predict observable properties and the evolution of QSOs 
at z > 6 (§4). In section 5, we discuss how the distant 
QSO J1148+5251 at z = 6.42 forms and predict the na- 
ture of the early phase of this QSO. Section 6 is devoted to 
our summary. Unless otherwise stated, all results shown 
below refer to the currently favored A cold dark matter 
model with SIm = 0.24, Ha = 0.76, h = 0.73, fib = 0.042, 
0-8 = 0.74 and n = 0.95 (Spergel et al. 2007). 

2. MODELS 



We here briefly describe a coevolution model of SMBHs 
and CNDs (KW08). We assume that dusty gas is supplied 
to a region around a central SMBH at a constant rate of 
Msup from a host galaxy whose surface density Shost, in- 
cluding the gas and stellar components (Fig. 1 in KW08). 
The accumulated gas forms a clumpy CND around a cen- 
tral SMBH with Mbh, which is vertically supported by 
turbulent pressure via SN explosions (WN02; Vollmer & 
Beckert 2003; Vollmer ct al. 2008; Collin & Zahn 2008; 
Wada et al. 2009). We here assume the isothermal cold 
gas dominates mass (Tg = 50 — 100 K) in the CND since 
the molecular and dust cooling is effective (e.g., Wada & 
Tomisaka 2005; Wada et al. 2009). 

2.1. Turbulent pressure-supported CND 

On the vertical structure of the CND, we assume that 
the turbulent pressure associated with SN explosions is 
balanced to gravity, g caused by 

Pg{r)vt{r) = Pg{r)g{r)h{r), (1) 
where Pg{r), Vt{r) and h{r) are the gas density, the tur- 
bulent velocity and the scale height of the disk, respec- 
tively. Here, the gravity, g(r) is obtained as g{r) = 
GAlBuh/r^ + 7rG(i;disk(r) + Shost) where Sdisk(r) is the 
surface density of baryonic components (the gaseous mat- 
ter and stars). The geometrical thickness is determined by 
the balance between the turbulent energy dissipation and 
the energy input from SN explosions as follows. 



riSiS!S^{r)EsN, 



(2) 



tdis{r) h{r) 
where the dissipation timescale of the turbulence, tdis{r) = 
h{r)/vt{r), EsN is the total energy (10^^ erg) injected by 
an SN and ?7sn is heating efficiency per unit mass which 
denotes how much energy from SNe is converted to kinetic 
energy of the matter. Recent observational studies on the 
relationship between star formation rate and gas surface 
density in nearby galaxies suggest that the star formation 
rate is approximately proportional to the gas density for 
high density region (nn > lO^cm"^), i.e., = C^p^ with 
n ~ 1 (e.g., Bigiel et al. (2008)). Here C* is the star for- 
mation efficiency. Theoretical studies based on numerical 
simulations of the inter stellar medium (ISM) also sup- 
port this proportionality for the high-density end (Wada 
& Norman 2007; Dobbs & Pringle 2009; Krumholz et al. 
2009). Since we are interested in the evolution of dense 
CNDs with riH > 10^ cm^'^, we here suppose 5* = C*/9g. 

The star formation time scale or C^^ is of the order of 
an orbital period in a quasi-steady, self-regulated galac- 
tic disk, in which the Toomre Q-value is roughly unity. 
However, observations suggest that the local star forma- 
tion rate for a given gas density varies widely distributed 
by one to two orders of magnitude (Komugi et al. 2005; 
Bigiel et al. 2008, see also Fig. 5 and discussion in §3.1). 
Based on a theoretical model of the ISM and star for- 
mation, Wada & Norman (2007) claimed that this wide 
variety could be caused by a wide range of star forma- 
tion efficiency. In general, the star formation rate (or ef- 
ficiency) on small scale is not simply determined only by 
gas density (see e.g. Kawamura et al. 2009 for molecular 
clouds in LMC). The physical background of this variety 
is still unclear even in our Galaxy. This is also the case 
in the galactic central region. Although many local spiral 
galaxies show molecular gas concentration in the central 
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sub-kpc region (e.g. Komugi et al. 2008), the star form- 
ing activity is not uniquely determined by the gas density. 
In fact, the local star formation rate is correlated with 
a fraction of high density molecular gas (Muraoka et al. 
2009). Moreover, at high rcdshift, the star formation cfh- 
ciency in the galactic central region is hardly determined 
observationally and theoretically. The circumnuclear gas 
disk is not necessarily in a self-regulated state, since the 
rotational time scale is comparable to the life time of mas- 
sive stars or smaller. Therefore under frequent mergca's of 
small gas-rich galaxies at high redshift, which would be a 
main cause of mass supply toward a galactic center, the 
circum-nuclear region might be far from in an equilibrium. 
Considering all these uncertainties, we here assume that 
the star formation time scale is a free parameter, and ex- 
amine the results for a wide range of C* . 

By using equations (1) and (2), Vt{r) and h{r) can be 
obtained (see also eqs. (4) and (5) in KW08). For given 

1 /2 1 /2 

r and Mbh, it is found that Vt <xCJ and h (xCj (see 
also WN02). The velocity dispersion of molecular hydro- 
gen in CNDs is in fact positively correlated with the star 
formation rate (Hicks et al. 2009). 

2.2. Two modes of gas accretion in CND 

We suppose a kinetic viscosity as a source of angular mo- 
mentum transfer in the gas disk. Then, the mass accretion 
rate in a viscous accretion disk is given by 

M(r) = 2^j.t(r)Sg(r) ^^^^^^ 



dim 



(3) 



where the viscous parameter is I'ti'r) = avt{r)h(r), Sg 
is the surface density of the gas component in the CND 
and the angular velocity il(r) is given by the radial cen- 
trifugal balance. Here a{< 1) is a constant in time. The 
inner radius of the CND (rin) is determined by the dust 
sublimation radius, i.e., = 3pc (Mbh/IO^M©)^/^ (see 
KW08 in details). At rin we assume the CND connects 
with the steady accretion disk. With respect to the sta- 
bility, we adopt Toomre's stability criterion, i,e., when the 
surface density of gas in the CND, Eg is higher (lower) 
than the critical surface density, Ecrit the CND is grav- 
itationally unstable (stable). The critical surface den- 
sity is obtained as Scrit('') = t(^)cs/7i'G, where K(r) = 
4f2^(r) -|- 2Q{r)rdQ{r) / dr is the epicyclic frequency and Cg 
is the sound velocity. The critical radius Tc is determined 
by the Toomre criterion, that is. Eg (re) = Ecrit(''c)- In this 
picture, we naturally expect two modes of gas accretion 
rate as follows: m,ode (i): If < Tin (i.e., the CND is fully 
gravitationally unstable), then the CND is geometrically 
thick due to SN energy feedback, and as a result we have a 
large accretion rate. We here suppose asN = 1 motivated 
by numerical simulations demonstrated by WN02. mode 
(ii): If fc > Tin, the scale height of the inner region would 
be much smaller than mode (i), because the scale height is 
determined by the thermal pressure, Pg{r) = pg{r)gh{r), 
where -Pg(r) = pg{r)c'^. Here Cg = {5kTg/'3mp}^^^ , where 
k and nip are the Boltzman constant and the proton mass. 
In mode (ii), the magneto-rotational instability could be a 
source of turbulence, but the turbulent velocity is compa- 
rable to or even smaller than the sound speed (e.g., Bal- 
bus & Hawley 1991; Machida et al. 2000). As a result, 
the accretion is less efficient than the mode (i). We as- 
sume aiviRi = 0.5 and t;t = Cg in mode (ii). Note that the 



adopted q;mri = 0.5 is relatively larger than that derived 
by numerical simulations (e.g., Machida et al. 2000), but 
we confirmed that the following results do not change at 
all, even if we use smaller ttMRi (e.g., ckmri = 0.01). 

2.3. CoevoluUon of SMBHs and CNDs 

Our main purpose is to evaluate growth of SMBH, the 
star formation rate and gas mass in the CND. We here 
focus on the time dependence of characteristic radius in 
the disk (rc, Tin, and Tout), instead of solving the evolu- 
tion of radial structure of CND. The radial surface den- 
sity distribution of CND is assumed to be Edisk(?') = 
Sdisk,o(7'/?'out)~''' where Tout is the outer radius of the 
CND. Hereafter, we assume 7=1 although the depen- 
dence of 7 on Mbh is weak for < 7 < 2 (see §3.2 in 
KW08). In this model the gas supplied from the host 
galaxy is eventually consumed to form the SMBH and 
stars. Thus, the time-evolution of the gas mass in the disk, 
Mg = J^*^""* 27rr' Eg (r')dr', is simply given by the mass con- 
servation: 



Mg(t) 



[M,up(t') - M,{t') - MBuit')]dt', (4) 



where Msup{t), M^{t) and MBuit) are the mass-supply 
rate from hosts, the star formation rate, and the growth 

rate of SMBH, respectively. Here we ignore the mass loss 
from stars and from CNDs due to the starburst wind. 
The time-evolution of SMBH mass MBH(i) is obtained as 

Mbh(^) = A/"BH,sccd + /q MBH{t')dt' , where we assume the 
mass of seed BHs, AfsH.sccd = 1O^M0, as end-products of 
the first generation stars (e.g.. Fryer et al. 2001; Heger et 
al 2003; Omukai & Palla 2003; Yoshida et al. 2006). For 
M 



sup J 



we can assume any function for the mass supply 
rate, but here we simply take a step function as the first 
attempt as Msup{t) = const for t < ts^p, while Msup(i) = 
for t > tsup where tgup is a period of the mass-supply from 
host galaxies. Thus, key parameters of the SMBH growth 
and the state of CNDs are M^^p , ^sup and . 

The growth rate of SMBH, i.e., Mbh, is not necessar- 
ily equal to the mass accretion rate at the inner bound- 
ary, M(rin). The growth rate could be limited by the 
Eddington accretion rate, MEdd where MEdd = L^dd/c^ 
and I/Edd = 47rcG'MBH'7ip/o"T with ctt being the Thom- 
son cross section. However, based on radiation hydrody- 
namic simulations, it is also claimed that Mbh > MEdd) 
i.e. super-Eddington accretion, could be possible since 
both the radiation field and mass accretion flow around a 
SMBH is non-spherical (Ohusga et al. 2005, 2007; Ohsuga 
& Mineshige 2007). Thus we here consider two possible 
models for Mbh- 

(i) Mbh = M(n„). (5) 
The BH grows with a super Eddington rate if M(rin) > 

MEdd- 

(ii) Mbh = ein^Edd, (6) 
where cbh = 0.1 is the energy conversion efficiency. We 
here assume that there is a mass "outflow" with a rate, 
-^outflow = -M(rin) — Mbh- Hereafter we call the model (i) 
and model (ii) as the super-Eddington growth model and 
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the Eddington-limited growth one, respectively. We first 
show results of the super-Eddington accretion in §§3.1 and 
4.1. Results based on the Eddington-limited accretion are 
discussed on §§3.2 and 4.2. 

Finally, we should stress what happens if we use the 
non-hnear Kennicutt -Schmidt's law, i.e., 5* = C^tg^pl-^, 
where tg is the free-fall time scale. It is found that the 
BH growth rate is then described by Mbh oc C^pg oc Mg 
because of the turbulent viscosity ft oc C^, while the star 

3/2 

formation rate scales as Af* oc Mg . Thus, the gas sup- 
plied from host galaxies turns into stars in the CND more 
efficiently, compared with the case for 5* oc pg (see eq. 
(4)). It implies that the final BH mass becomes smaller 
due to the efficient star formation, if we adopt non-linear 
Kennicutt-Schmidt's law. 
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The final SMBH mass, Mbh, final against the diff'er- 



tsup = 10 yr, 10 yr and 10 yr. The red lines are the results of the 
super-Eddington growth model, while the black lines are ones of the 
Eddington-limitted growth model. The blue line denotes that all the 
supplied mass is used to grow the SMBH, i.e., Afgjj (j^a^i = Mgup. 
Here we assume C» = 3 X 10~* yr~^. 



3. NECESSARY CONDITIONS OF FORMATION OF SMBHS 
AT Z > 6 

Based on the coevolution model of CNDs and SMBHs 
(§2 and KW08), we examine conditions of the formation 
of SMBHs with Mbh > IO^Mq at z > 6. We first show 
results of the super-Eddington growth model in §3.1, and 
those of the Eddington-limited growth model in §3.2. 

3.1. Super-Eddington growth model 

Figure 1 shows the final BH mass, A/bh. final as a func- 
tion of the total mass accreted from hosts, Mgup for the 



different period of mass supply from hosts, tgup assuming 
C* 3 X 10^* yr^^. This value is comparable to the star 
formation efficiency of high- z galaxies (e.g., Tacconi et al. 
2006; see also Fig. 5). The effect of C* will be discussed be- 
low. The red lines are the results of the super-Eddington 
growth model for tsup = 10^ yr, 10® yr, and 10^ yr. The 
blue line denotes that all the supplied mass is used to 



grow the SMBH, i.e., A/bh. final = Mgup- It is found that 
ATbh, final ~ (0.01 — 0.1)A/sup- This low efficiency is in- 
dependent of the period of the mass-supply from hosts, 
isup and the mass supply rate of hosts, Mgup. In order to 
form SMBHs with > lO^M©, plenty of infalling gas from 
host galaxies, i.e.. A/gup > IQ-^^M© is requested. We con- 
sider why the mass ratio, il^BH.finai/Afsup is too small as 
we showed. From the BH growth rate for the high ac- 
cretion rate (§2.2 and see also eq. (12) in KW08), the 
accretion energy onto a central BH, i.e., O.SMeHi^in where 
Win = \/ G Mbh /rim can be obtained as 



^AfBH^f, 



= 2/^asN?7sN-BsNM* 



(7) 



where = [3(14 - 47)(8 - 47)/16] (nn/rout)"-^- Thus, 
equation (7) indicates that asN and t^sn are important pa- 
rameters to determine the conversion efficiency from the 
energy input from SN explosions (« EsnM*) to the ac- 
cretion energy onto a central BH (« MBH^'fn)- Using the 
fiducial values of 7 = 1, the ratio of the SMBH growth 
and star formation rate of CNDs is given by 



Af. 



= 0.3q?sn 



10-3Af 







VlO^AfG 



(8) 



This indicates that larger Mgup (or Afdisk) leads to smaller 
A^BH.finai/Afsup- This cau be understood as follows. As 
Afsup becomes larger Mbh increases. The larger BH mass 
causes a smaller scale height and a smaller turbulent ve- 
locity because the CND is vertically supported by the tur- 
bulent pressure (eq.(l)). Thus, the growth rate of BHs 
(oc vt oc vth) decreases as Mgup (or Afdisk) is larger. Note 
that A/disk = Afg -I- M* is comparable to A/gup. If we focus 
on the case for Mgup > IO^Mq, it turns out that the BH 
growth is always smaller than the star formation rate in 
the CND, i.e., Mbh < 0.3M,,, assuming that asN = 1 and 
'/SN = 10"'^ Mg^ 0. In other words, the star formation 
rate always overcomes the mass accretion rate onto BHs 
in this model. From the eq. (8), at given Mgup (or A/disk) 
the ratio A/bh/A/* is basically determined by two quan- 
tities, i.e., asN and ?7sn which control turbulent viscosity 
and efficiency of SN heating, respectively. 
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* The numerical simulations have shown ag^ 1 (e.g., WN02). On the other hand, the parameter ?7sn is expressed as t^sn = egf,j/gN, where 
egt-j and /gjvj are the efficiency with which SN energy is transferred to the gas in the CND, and the number density of SNe per solar mass of 
the star formation, respectively. In this paper (also in KW08), we assume r;sN = lO^^Mg"*^, that is, egN = 0.1 (e.g., Thornton et al. 1998; 
Wada & Norman 2002; Wada et al. 2009), and /gN = 10~^Mq^ which is expected for Salpeter initial mass function (IMF) with low-mass 
cutoff being mi = O.IMq. Note that the heating efficiency jygjsj must be larger than 10"'^ M^^^ in order to be satisfied with Mbh > M». 
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Fig. 2. — Time evolution of the mass of BH for different C*. 
The right blue, magenta, green, blue and red lines correspond to 
C* = 3 X lO-'^yr-l, 3 x 10-''yr-l, 3 x 10"* yr"!, 3 x IQ-^yr-i 
and 3 X 10"^'' yr~^, respectively. Here we assume Maup = IO^'^Mq 
and tsup = 10^ yr. The shaded region represents the allowed region 
for the formation of SMBHs with Mbh > Mq at 2 > 6. The 
dot-dashed vertical line corresponds to 2 = 6 if the redshift Zi at 
which the accretion onto a seed BH starts is = 25. 
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Fig. 3. — Same as Fig. 2, but tsup = 10* yr. The lines are same 
as Fig. 2. 

Next, we examine the dependence of the star forma- 
tion efficiency, C* = M*/Mg on the mass of SMBHs, 
MBH(i)- Figure 2 shows the evolution of the SMBH 
mass, Mbh(^) for the different star formation efficiency, 
C*. We here assume the total mass accreted from hosts, 
Msup = IO^^Mq and a period of the mass-supply from 
hosts, tsup = 10^ yr. The redshift Zi when the accretion 
onto a seed BH got started is supposed to be Zi = 25. 
Thus, the growth timescale of SMBHs must be less than 
8 X 10^ yr (vertical dot-dashed hne) to form the SMBH 
at z > 6. The allowed region to form a SMBH with 
-^^BH > IO^Mq at z > 6 is shown as shaded. We find 
that MBH(i) significantly depends on the star formation 
efRciency, C* . If the star formation efficiency is extremely 
low (e.g., < 3 X lO-^r-i), the timescale of SMBH 
growth is larger than the Hubble time at z = because of 

^growth ~ 

/vth - (AC*)-! (ggg (;l4) KW08). 
Thus, A/rh, final cannot be larger than ~ lO^M© until 
z — 6 (the red line in Fig. 2). For the extremely high 
efficiency (e.g., C* > 3 x 10"^ yr"-'^; light blue line in 
Fig. 2), on the other hand, the gas supplied from the host 
is mostly consumed by the star formation (Msup < M^), 
thus the final SMBH mass is smaller than IO^Mq. There- 
fore, in order to form a SMBH with A/bh > W^Mq 
at z > 6, relatively high star formation efficiency, i.e., 
10~^yi~^ <Cf:< IQ-^yr--'^ is required. 

We show the evolution of A/bh(0 for the case of isup = 
IQSyr (Fig. 3) and t^up = 5 x lO^yr (Fig. 4) respec- 
tively, to elucidate the dependence of a period of the 
mass-supply from hosts <sup on M-QY{{i)- Comparing Fig. 
2 with Fig. 3, the final BH masses for igup — 10* yr 
are almost the same as that for igup = 10^ yr (Fig. 2), 
but the behavior of A/BH(i) for isup = 10* yr is different 
from that for tsnp = 10^ yr at given C* = 3 x yr"^ 
(the hght blue line). This is because for ts^p = 10* yr 
the mass supplied rate is higher than the star formation 



rate (Afgup > M^,). Thus, the final BH mass can be as 
large as ~ lO^Af©. As a result, the allowed region is 
10-*yr-i < C* < 10-6 yr-i for U^p = 10* yr. The upper 
value C* is determined by the local gravitational free-fall 
time scale, ts « (Gridisk^p)^^^^ ^ 10^ yr for the averaged 
density of CNDs, i.e., ndisk = lO'^^^cm^'^, because of the 
star formation time scale = C^^ > tg. We also found 
that the required range of for tgup = 10^ yr is the same 
as those for tgup = 10* yr. Figure 4 shows that there is no 
solution to form the SMBH with AfsH > W^Mq at z > 6 
for tsup = 5 X 10^ yr because of Afgup < A/* for the wide 
range of C*, i.e., 3 x 10'^° yi'^ < < 3 x IQ-^yr-^. 
In other words, this indicates the rapid mass supply from 
the host to the central tens pc region (Afsup > 20Afo yr--^, 
Msup = lO^Afg) is also necessary to form the SMBHs at 
z > 6. 
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Fig. 4. — Same as Fig. 2, but for tsup = 5 X 10^ yr. The lines are 
same as Fig. 2. 

Figure 5 shows the SFR as a function of the surface 
density of gas in CNDs, Sg. We estimate the surface gas 
density of CNDs as Sg = M^j-Kr^^^. Here we take Afg and 
Tout from our calculations. The size of the CND increases 
with time but finally Tout is ~ 2 x 10'^ pc. For the SFR and 
Sg in the shaded region, SMBHs with AfBH,finai > lO^Af© 
at z > 6 can be formed as a result of considering the de- 
pendences of C* and Eg. The observational data of nearby 
normal galaxies (Komugi et al. 2005) and nearby star- 
burst galaxies (Kennicutt 1998) are plotted in Fig. 5. The 
two plots represent the average of massive star forming 
submillimeter galaxies (SMGs) at z ~ 2.5 with SFR w 
lO-^Afgyr-ipc-^ and Sg « 5 x IO^ATqpc-^ (e.g., Tac- 
coni et al. 2006) and a QSO at z = 6.42 (SDSS J1148) with 
SFR « lO-3M0yr-ipc-2 and Sg « IO^'ATq pc-^ (Walter 
et al. 2009). The required C, for high-z SMBH forma- 
tion is comparable to SMGs and QSO hosts a,t z — 6.42, 
although this is higher than that of low-z normal galaxies 
and starburst ones. Thus, the predicted high star forma- 
tion efficiency, i.e., 10-* yr-^ < C* < IQ-^yr-^ actually 
occurred in the high-z universe. 
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Fig. 5. — Comparison between the observed surface star formation 
rate (SFR) in terms of the surface gas density (Sg) and the allowed 
region for the formation of SMBHs with Mbh > 10^ Mq at 2 > 6 
(the shaded region). The blue lines represent theoretical predic- 
tions for different C«, assuming Msup = 10^^ Mq and tsup = 10* yr. 
Black filled dots are starburst galaxies in Kennicutt (1998) and the 
black open dots are normal galaxies in Komugi et al. (2005). Green 
point denotes the average dusty submillimeter galaxies (SMGs) at 
2 2.5 (Tacconi et al. 2006) and red point represents ~ kpc scale 
starforming region around the z = 6.42 SDSS J1148+5251 (Walter 
et al. 2009). 

In summary, in order to form the high-z SMBH (i.e., 
z > 6 and Mbh, final > IO^Mq) by super-Eddington 
gas accretion we need; (i) abundant mass supply from 
the galactic scale (~ kpc), i.e., Mgup > IO^^Mq, (ii) 
rapid mass inflow from the galactic scale {'^ kpc), i.e., 
tsup < 5 X 10^ yr and (iii) high star formation efficiency, 
i.e., fO-^yr-i < < 10-^yr-\ 

3.2. Eddington-limited growth model 

We here discuss the results of Eddington-limited growth 
model (the black lines in Fig. I). Figure 1 shows that 
the Eddington-limited growth model leads to much lower 
efficiency for i^up = lO'^yr and 10^ yr, i.e., Mbh. final ~ 
(10"'' - 10-2)Msup. But, the efficiency for t^up = lO^r is 
comparable to that of the super-Eddington growth model, 
i.e, Mbh, final ~ (O.OI — O.I)Msup. These results can 
be understood in terms of the ratio of the mass-supply 
time scale, igup and the Eddington time scale, t^dd = 
eBHO'Tc/(47rGmp). The black hole grows exponentially at 
the late time {t ^ <Edd ~ 5 x f 0'' yr for ebh = 0.1) because 
of Mbh ck e*/*'^''''. If tgup is comparable and/or shorter 
than iEdd (e.g., isup = 10 yr and 10* yr), the mass-supply 
from a host galaxy terminates before a main growth phase 
of BHs. Thus, the BH growth is greatly suppressed by the 
AGN outflow (see eq. (6)), and then the final BH mass can 
be 1-2 orders of magnitude smaller than that of the super- 
Eddington growth model as seen in Fig. 1. On the other 
hand, the effect of AGN outflows in high accretion phase 
is not significant for igup = 10^ yr because of tgup 3> ^Edd- 
Thus, the final BH mass is comparable to that of the super- 
Eddington growth model. In order to form a SMBH with 
> 1O®M0, it is found that the super-Eddington growth is 
needed when tgup is shorter than ~ 10^ yr. Note that the 
dependence of C* and igup on the final BH mass are the 
same as those for the super-Eddington growth model. 

4. OBSERVABLE PROPERTIES AND EVOLUTION OF QSOS 

AT Z > 6 



Based on the present model in §3, we here discuss ob- 
servable properties and the evolution of QSOs at z > 6, 
in order to compare them with the current and future ob- 
servations. It is still unclear when the accretion onto a 
seed BH got started in the high -z universe. Recently, 
Johnson & Bromn (2007) mentioned that the radiation 
feedback from the first stars may deplete the gas in the 
central region of the proto-galaxy and may delay the BH 
accretion by up to ~ 10* yr. Thus, we explore the rela- 
tionship between the evolution of QSOs at z > 6 and the 
redshift Zj when the seed BH begins to grow. To this end, 
we examine the evolution and formation of high-z QSOs 
at z > 6 for the two scenarios of high-z QSO formation; 
(a) the early growth of SMBHs, i.e., Zj = 25. This cor- 
responds to the supply period tgup = 7 x 10* yr and (b) 
the late growth of SMBHs, i.e., Zi = 8. This corresponds 
to the supply period tgup = 10* yr. We here use the fidu- 
cial values, i.e., M^up = 10"Mq and C* = 3 x 10-*yr-\ 
for which a SMBH with Mbh > 10^ Mq at z > 6 can be 
formed. We define the QSO phase as tgup < t < tqso 
where tqso is the time when the AGN luminosity equals 
to the threshold luminosity of QSOs, Lgscth- We assume 
-^Qso.th = lO^^ergs"^. The proto-QSO phase is defined 
as the early phase of a growing BH {0 < t < isup)- 
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Fig. 6. — Time evolution of optical depth of CNDs for UV-band 
in the edge-on view (t). The blue line corresponds to the case for 
2i = 25 (or 4sup(a) = 7 X 10® yr). The red line denotes the case for 
2i = 8 (or tsup(b) = lOVr)- 

Figure 6 shows the evolution of optical depth of CNDs 
for UV-band in the edge-on view, r(t) defined as T{t) = 
Ir.y XdPg{r)dr ~ XdMg{r)/[2TTrl^^{r)hout{r)] where Xd 
is the mass extinction due to the dust and hout{f) is 
the scale height at rout(?')- The mass extinction coef- 
ficient Xd is given by Xd = nd<Jd/Pg with the number 
density rid and cross section ad- In this paper, we as- 
sume Xd = 103cm2g-i(ad/0.1/xm)-i(ps/gcm-3)(Z/2'Q), 
where ad is the grain radius, ps is the density of solid ma- 
terial density within the grain (e.g., Spitzer 1978; their 
§9.3), and Z is the metallicity of gas, which are fixed at 
Od = 0.1/im, Ps = lgcm~^, and Z = Zq. We find that 
CNDs are always optically thick during the evolutionary 
path of bright QSOs for both cases (i.e., Zj = 8 and 25). 
The present model predicts the existence of an optically 
thick CND even in a proto-QSO phase. At late times 
{t > 10^ yr), the properties of CNDs (e.g., Mg, Tout, and 
/lout) are independent of time because of no fueling and 
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no star formation in this phase. Thus, the optical depth 
is constant in time. 

In order to investigate the evolution of the AGN lu- 
minosity and the nuclear-starburst luminosity, we define 
these luminosities as follows: Following Watarai et al. 
(2000), the AGN luminosity can be given as a function 

of mBH = MBn/MEdd 



z=25.0 



LAGNit) 



2 1 + In ■ 



tBH(t) 

20 



^Edd 



(t) ;mBH(t) > 20, 



(^^)LEddW ;mBH(t) 



< 20, 



(9) 

where ttibh = Mbh/ MEdd- On the other hand, the nuclear 
starburst luminosity isB(i) can be obtained as 



LsB{t) = 0.14e*M*(t)c2 



(10) 



where = 0.007 which is the energy conversion efficiency 
of nuclear fusion from hydrogen to helium. The results 
of the super-Eddington growth model are shown in §4.1. 
Results based on the Eddington-limited accretion are dis- 
cussed on S4.2. 



tsup for Zi = 8, while it reaches the maximum much 



10 < t,np, for 



4.1. Super-Eddington growth model 

Figure 7 (a) and (b) shows the evolution of bolometric 
luminosity of AGN, Lagn (red lines), nuclear starburst 
luminosity, LsB (green lines) and Eddington luminosity, 
^Edd (blue lines) for two scenarios of high-z QSO forma- 
tion. The shaded region shows QSO phase. We find that 
the evolution of QSO luminosity depends on the redshift 
Zj at which accretion onto a seed BH is initiated, although 
the final BH mass, A^bh, final and the evolution of BH mass, 
Mbh arc independent of Zj (or tsup) at given C* (see Fig.l, 
2 and 3). The time when Xagn becomes maximum is at 
t 

earlier than the QSO phase, i.e., t 
Zj = 25. In addition, the maximum luminosity of AGNs 
for Zj = 8 is an order of magnitude larger than that for 
Zi = 25. These results indicate that the evolution of the 
QSO luminosity is sensitive to the redshift Zj when the 
accretion onto a seed BH is initiated. In other words, it 
is easier to build-up bright QSOs if the gas is supplied 
rather late on a relatively short timescale. This can be 
understood as follows. For Zj = 8, the mass supply time 
scale, tsnp is shorter than the star formation time scale, 
oc C^^. Thus, the gas in the CNDs can accrete onto a 
BH more efficiently, compared with the case for Zi = 25. 
On the contrary, if the average gas supply rate is smaller in 
a longer period, (i.e. Zi = 25), the supphed gas is mainly 
consumed to form stars not to form a SMBH because of 
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Fig. 7. — (a) Time evolution of the AGN luminosity, LAGN(i)i 
that of the nuclear starburst luminosity, Lggit) and that of Ed- 
dington luminosity, -LEdd(t)i assuming 

3 X 10-«yr-i and = 25 (or isup = 7 X lO^yr). We define 
QSO phase as tsup < t < tqso (the shaded region), where tqso 
is the time when the AGN luminosity equals the threshold luminos- 
ity of AGN, iqscth = Iff'^erg/s. For A^bh growth, we assume 
the super-Eddington growth model, (b) Same as (a), but for z; = 8 
(or tsup = 10* yr). 

We also find that the AGN luminosity is sub-Eddington 
with LAGN/^Edd ~ 0.1 — 1 in the QSO phase, while it is 
super-Eddington with Lagn/ LEdd « 1 — 10 in the proto- 
QSO phase (see Fig. 7(a) and (b)). On the other hand, the 
luminosity for the QSO phase is comparable to than that 
in the proto-QSO phase. This indicates that the Edding- 
ton ratio, -ZjAGN/iEdd, varies two orders of magnitudes for 
a given AGN luminosity. This implies that the BH mass 
derived from only AGN luminosity leads to a large uncer- 
tainty. 

The nuclear starburst (SB) luminosity for the proto- 
QSO phase is larger than that for the latter half of QSO 
phase, because the SB luminosity increases with time in 
the proto-QSO, while it decreases with time in the QSO 
phase. This trend is prominent for Zj = 8 rather than for 
Zi = 25 (see Fig. 7 (a) and (b)). In this model, the enhanced 
star formation in a BH growing phase (i.e., a proto-QSO 
phase) is expected. At low-z universe, Sani et al. (2009) 
discovered that the star formation activity in narrow line 
Seyfert galaxies (NLSls) is larger than that in broad line 
Seyfert galaxies (BLSls) of the same AGN luminosity. If 
NLSls are the early phase of BLSIs, these observations 
are consistent with our predictions. We also find that the 
nuclear SB luminosity is comparable to the AGN luminos- 
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ity in both the proto-QSO and QSO phases. Since young 
stars are buried in the optically thick CNDs (see Fig. 6), 
the intense nuclear starburst with Lsb ~ lO'*^ — 10'*^ erg/s 
can be observed in the infrared band. 

The lifetime of QSO is tqso ~ 1 x 10* yr for zi = 25 and 
T'QSO ~ 3 X 10® yr for Zi = 8. These values of tqso are con- 
sistent with the previous studies (e.g., Martini et al. 2003a, 
b; Kawakatu, Umemura & Mori 2003; Granato et al. 2004; 
Shen et al. 2007; Li et al. 2008). From comparison be- 
tween two models, we find that the QSO phase is longer in 
more luminous QSOs. We should note that tqso is larger 
for less luminous QSOs, e.g., tqso « 2 x 10® yr (case (a)) 
and Tqso ~ 5x10® yr (case (b)) for iqscth = lO^^ergs"^ 
This has already been pointed out by Hopkins et al. (2005) 
and Li et al. (2007). 

Figure 8 shows the time evolution of the BH mass, 
M-Q}i{t), the gas mass in the CND, Mg{t), the stellar 
mass in the disk, {t) and the total supplied mass from 
hosts. The evolution of Mg(i) is closely related to that of 
LAGN{t) because the CNDs become gravitationally sta- 
ble as Mg decreases (see §2). In the QSO phase, we 



find that Mn- 



and Af^ 
(b)) 



IOOMbh 



(3 



- 5)Mbh = 5 x 10^ - 1 X 101°Mq, 
10"Mq for Zi = 25 (Fig. 8 (a)), 
(3 - 30)AfBH = 5 X 10^ - 5 X IO^OM© and 
(10 - 100)Mbh = 10^° - 10"Mo for = 8 (Fig. 8 
Thus, it is predicted that there exists a stellar rich 
massive CND around SMBHs, i.e., Mg w O.lAfdyn where 
Afdyn is the dynamical mass of the CND plus BH system. 
On the other hand, in the proto-QSO phase it is found that 
Afg/MBH = 10- 10^ for Zi = 25 and Afg/A^BH = 10^ - 10^ 
for Zi — 8. Thus, we predict that the proto-QSOs have 
a gas rich CND around SMBHs, i.e., A^g ~ Afdyn- At 
late times {t > 10^ yr) , a very massive stellar nuclear disk 
(Af* ~ IQ^^Mq) remains as a by-product of the formation 
of SMBHs with A^bh > 10^ Af©. Since the size of torus 
(''out) grows up ~ 2 kpc in this phase, the circular velocity 
of stellar components in the CND is ~ 300kms~^. This 
is comparable to the velocity dispersion of giant elliptical 
galaxies (e.g., Tremaine et al. 2002). Comparing the case 
of Zi = 25 with that of Zj = 8, the larger amount of gas 
in the CND for Zi — 8 exists just before the QSO phase 
(see Fig. 8 (a) and (b)). As a result, the higher mass 
accretion rate (A/bh) is achieved, and the brighter AGN 
luminosity is sustained for a longer timescale. This is why 
the QSO lifetime for Zj = 8 is longer (see Fig. 7 (a) and 
(b)) than that for Zi ~ 25. In addition, for Zi = 25 the 
stellar mass Af* becomes larger than the gas mass Afg at 
^max ~ 10® yr. Thus, the AGN luminosity does not in- 
crease after tj^ax ~ 10® yr (see Fig. 7 (a)), even when the 
gas continues to infall from hosts into CNDs. This sup- 
plied gas is mainly consumed to form stars, not to form 
SMBH. 

Based on the present results, we can predict the observ- 
able properties of proto-QSOs as follows; (i) a gas rich 
CND with the mass ratio Afg/AfBH > 10^"^, which is 
larger than that of QSOs, i.e., Afg/A/BH < 10. (ii) the nu- 
clear SB luminosity is slightly higher than that of QSOs. 
Thus, in order to explore proto-QSOs, it is essential to 
estimate the gas and stellar mass in CNDs, and the nu- 
clear SB activity. Jiang et al. (2008, 2009) have recently 
discovered two magnitudes fainter QSOs at z 6 than 
bright QSOs discovered in previous survey (e.g.. Fan et al. 



2003). Moreover, Kurk et al. (2009) found that two faint 
QSOs have very narrow broad emission lines (Mg II and 
CIV), which may imply small BH mass and high Edding- 
ton luminosity ratios. Thus, these faint QSOs may be good 
candidate of proto-QSOs with a super-Eddington mass ac- 
cretion flow. It is interesting to examine the presence of 
gas rich CNDs for these high-z faint QSOs using ALMA. 
But, we should mention that it is hard to judge whether 
the proto-QSOs are super-Eddington objects by only these 
observable properties of CNDs. Thus, we will really need 
to explore the direct evidence of super-Eddington objects. 
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Fig. 8. — (a) Time evolution of the BH mass, AfBH(<), the 
gas mass in the CND, Afg(t), the stellar mass in the disk, 
and the total mass supplied from hosts, Afsup(t) for Zi = 25 (or 
^sup = 7 X 10* yr). The shaded region denotes the QSO phase. For 
A^BH growth, we assume the super-Eddington growth model, (b) 
Same as (a), but for Zi = 8 (or tsup = 10* yr). 



4.2. Eddington-limited growth model 

In §4.1, we have shown the evolution of luminosity 
and mass for the super-Eddington growth model (see Fig. 
7 and Fig. 8). In order to compare with the super- 
Eddington growth model, we examine these evolutions 
for the Eddington-limited growth models with Zi = 8 (or 
tsup = 10® yr) in Fig. 9, and Zi = 25 (or i^up = 7 x 10® yr) 
in Fig. 10. Figure 9 (top panel) shows that the peak 
AGN luminosity is ~ 10"*^ erg/s. This is one order magni- 
tude smaller than the threshold luminosity of QSOs (i.e., 
^QSO,th = 10*^ erg/s) and thus the QSO phase does not 
appear. This is because the BH growth is greatly sup- 
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pressed (see bottom panel in Fig. 9 and discussion on 
§3.2). Thus, the Eddington-limited growth model for 
Zi = 8 cannot explain the observed distant luminous 
QSOs. Also, the nuclear starburst activity dominates 
the AGN activity, and the peak nuclear starburst lumi- 
nosity is ^ 10'^''erg/s. Such an intense starburst would 
be observed in the infrared band because of the strong 
dust extinction. If this is the case, we may discover a 
hyper- luminous infrared galaxy without a bright AGN at 
z 7. Moreover, we predict a strong AGN outflow with 
-^outflow ~ lOM0yr~^ in the nuclear starburst dominated 
phase (bottom panel in Fig. 9). Here the total outflow 



mass is obtained by A^outflow = /q 



outflow 



{t')dt'. 
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Fig. 9. — (top panel): Same as Fig. 7 (a), but results for the 
Eddington-limitted growth model for Zi = 8 (or tsup = 10* yr). 
(bottom panel): Same as Fig. 8 (a), but results for the Eddington- 
limitted growth model for z; = 8 (or tnup = 10* yr). The total 
outflow mass (Moutflo„) is newly added. In this case, QSO phase 
(Lagn > iqscth) does not appear. 

Figure 10 (top panel) shows that the peak AGN lumi- 
nosity is comparable to ./jqso.th for the Eddington-limited 
growth model for Zj = 25. This is because the final 
BH mass is comparable to that for the super-Eddington 
growth model (see Fig. 1 and the discussion on §3.2). 
Thus, this model can marginally reproduce the observ- 
able properties of QSOs at z ~ 6. However, it is hard 
to form more luminous distant QSOs at z > 8 by the 
Eddington-limited growth as far as the radiative efficiency 
is eBH '^0.1. Thus, it is essential to discover the distant 
luminous QSOs by the future deep wide-field survey. Fig- 



ure 10 (top panel) also shows that the nuclear-starburst 
luminosity dominates the AGN luminosity in the proto- 
QSO phase (i.e., Lsb 3> Lagn)- This is different from 
the results of the super-Eddington growth model. On the 
other hand, the evolution of Mg {t) and M^, {t) are the same 
as the results of super-Eddington growth model (see Fig. 8 
(a)). Thus, it is predicted that a stellar rich massive CND 
exists in the QSO phase and for the proto-QSO phase there 
exists a gas rich massive CNDs (see bottom panel in Fig. 
10). 
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Fig. 10. — (top panel): Same as Fig. 7 (a), but results for the 
Eddington-limitted growth model for Zi = 25 (or tsup = 7 X 10* yr). 
(bottom panel): Same as Fig. 8 (a), but results for the Eddington- 
limitted growth model for Z[ = 25 (or tsup = 7 X 10* yr). The total 
outflow mass (A/outflow) is newly added. The shaded region denotes 
the QSO phase. 



5. APPLICATION TO THE DISTANT QSO SDSS 

J1148-I-5251 AT z = 6.42 

We here compare our predictions (§4) with a distant 
QSO, the wefl known and examined SDSS J1148-H5251 at 
z = 6.42. We summarize the properties of this QSO elu- 
cidated by using multi- wavelengths. The bolometric QSO 



luminosity is Lag^ 



10 erg s by near infrared ob- 



servations (e.g., Willott et al. 2003; Earth et al. 2003). 
The radio and sub-mm observations suggest vigorous star 
formation in ~ kpc with ^ 10'^ M0 yr~^ (Bertoldi et al. 
2003a; Carilh et al. 2004; Walter et al. 2009). In addition, 
the molecular gas in the central ~ 6 kpc is ~ 5 x IO^^Mq 
by using carbon monoxide (CO) emission (Walter et al. 
2003. 2004; Bertoldi et al. 2003b). However, we should 
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note that the observed emission hne might be sitting on 
top of a much broader hne, because the low signal to noise 
ratio of the current data is not enough to determine the 
emission with accuracy (Narayanann et al. 2008). 

The observed high luminosity with Lagn ^ 10'*^ erg , 
could correspond to the late growth of SMBH with Zi = 8 
(see Fig. 7 (b)) rather than the early growth of SMBHs 
with Zi = 25 (see Fig. 7 (a) and Fig. 10 (top panel)). If 
this is the case, the rapid and high mass accretion from ~ 
kpc scale (i.e., M^up « lO^Mgyr"^ and M^np « 10"Mq) 
is required. Thus, the presence of massive CNDs with 
Mg = 5 X 1O^°M0 and a luminous nuclear starburst in 
the infrared band with _Lsb ~ 10"*^ ergs" ""^ are predicted. 
These values are larger than the case of the early growth 
of SMBHs with Zi = 25, i.e., Mg < IQI^Mq and Lsb « 
2 X lO'^^ergs"^. The predicted CNDs can be observed 
through CO and hydrogen cyanide molecular emission us- 
ing ALMA even at z - 6 because of Mg « 10^"""Mq in 
the CND (e.g., Kawakatu et al. 2007). The PAH emission 
is a powerful indicator of the obscured starburst activity 
(e.g., Imanishi 2002). For 3.3 /im (the rest frame) PAH 
emission from CNDs at z = 6, ^ 3 x 10~^^Wm~^ line 
flux is expected at ^ 25 /im at the observed frame, assiim- 
ing that the ratio of 3.3 /xm PAH emission and the total 
infrared luminosity is w 10~^ for pure starburst galaxies 
(e.g., Mouri et al. 1990; Imanishi 2002). The Space In- 
frared Telescope for Cosmology and Astrophysics (SPICA) 
and the James Webb Space Telescope (JWST) through 
polycyclic aromatic hydrocarbon (PAH) emission could 
also detect dusty starbursts in CNDs of this QSOs be- 
cause the predicted value is above the line sensitivity limit 
(Swinyard et al. 2008). The present model can be verified 
by these future observations. 

On the basis of the present results, we will constrain 
the host mass of the z = 6.42 QSO SDSS J1148-I-5251. 
The mass accretion due to the tidal torque driven by the 
major and minor merger of galaxies is an essential process 
during the hierarchical formation of galaxies (e.g., Toomre 
& Toomre 1972; Mihos & Hernquist 1994, 1996; Saitoh & 
Wada 2004; Saitoh et al. 2008). Saitoh & Wada (2004) ex- 
amined the evolution of stellar and gaseous cores on a sub- 
kpc scale during the formation of a spiral galaxy by hierar- 
chical mergers, and found that the galactic core (< 300 pc) 
coevolves with the galactic dark matter (DM) halo of ^10 
kpc scale. The average mass ratio between the gas compo- 
nent in the core, Mcoro and the DM halo, Mdm is nearly 
constant, Mcore/-^DM ~ 0.04. If we assume that the galac- 
tic core corresponds to the CNDs, i.e., Mcore = -^sup; 
QSOs at z > 6 would be formed in a massive DM halo of 
Mdm > lO^^M© because of M^up > 10^° Mq. A massive 
halo of w 10^2" ^^Mq was also suggested by Li et al. (2007) 
to explain the low space density of bright QSOs at z 6 
(e.g.. Fan et al. 2001). Alternatively, Umemura (2001; 
hereafter UOl) proposed another mass accretion process, 
i.e., the mass accretion onto the galactic center from a 
galactic scale (~ kpc) via the radiation drag (see also 
Kawakatu & Umemura 2002). According to UOl, the max- 
imal mass accretion rate in the optically thick regime is 
given by Mdrag = ^^host/c^ « O.IMq yr-i(Lhost/lO^^L0), 
where Lhost is the luminosity of host galaxies. To re- 
produce the rapid mass infall, e.g., Mgup ~ 10M©yr^^, 
the hosts of proto-QSOs must have strong star forma- 



tion with ihost ~ 10^^X0. If this is the case, the stellar 
mass of the QSO host can be Mhost ~ 10^'^ Mq because 
of Mbh ~ O.OOlMhost (UOl; Kawakatu, Umemura and 
Mori 2003). These discussions indicate that the SDSS 
J1148+5251 form in massive dark matter halos and/or 
massive hosts (i.e., Mdm ~ 1O^^~-'^^M0 and/or Mhost ~ 
10^^ Mq), in order to build a SMBH with > lO^M© only 
by the gas accretion. This point should be tested by the 
future deep wide-field AGN survey. 

Moreover, we discuss the nature of the progenitor of 
SDSS J 1148+5251. As seen in Figure 7 (b), the AGN lumi- 
nosity is still bright with Xagn ~ lO^^ergs"^ at z ~ 7.5. 
It is hard to observe the AGN feature in the optical and 
X-ray observations because the AGN could be obscured by 
the dusty CND and dusty hosts (see Fig. 6). We predict 
that the progenitor can be super-Eddington objects (see 
Fig. 7 (b)) and observed as high-z ULIRGs such as SMGs 
with LiR w 10'*'^ ergs"-'^ because of Lagn ~ L^b- It is also 
predicted that the progenitor has a gas rich CND with 
the mass ratio Mg/Mdyn ~ 1- Such a gas rich CND with 
Mg ~ lOi°-"M0 can be detected by ALMA. 

6. SUMMARY 

Based on the physical model of supermassive black hole 
(SMBH) growth via gas accretion in ~100 pc scale circum- 
miclcar disks (CNDs) proposed by Kawakatu and Wada 
(2008), we investigate the formation of high-z: (z > 6) 
QSOs whose BH masses are Mbh > lO^M©. We show 
the necessary conditions to form QSOs at z > 6, in terms 
of the total gas mass accreted from host galaxies, Mg^p 
and star formation efficiency in the CND, C* = M*/Mg 
where and Mg are the star formation rate and the gas 
mass, respectively. Our main conclusions are summarized 
as follows. 

1. The required conditions for the formation of QSOs 

at z > 6 are as follows; (i) Rapid and large mass 
supply with Mgup > IO^'^Mq within 1 Gyr from 
hosts to CNDs, because the final BH mass is only 
1 — 10% of the total mass supplied from hosts. The 
fraction is determined by the efficiency with which 
SN energy is transferred to the cold gas, and the 
efficiency of angular momentum transfer due to tur- 
bulent viscosity, (ii) High star formation efficiency, 
i.e., 10~^yr~^ < C* < 10~^yr~^ is comparable to 
high-z starburst galaxies such as SMGs. (iii) Rapid 
mass inflow from hosts with a period of the mass- 
supply, tsup < 5 X 10^ yr. 

2. We find that if the BH growth is limited by the Ed- 
dington accretion, the final BH mass is greatly sup- 
pressed when the period of mass-supply from hosts, 
tsup, is shorter than the Eddington timescale. Thus, 
the super-Eddington growth is required for the QSO 
formation, while tgnp, which is determined by the 
efficiency of angular momentum transfer, is shorter 
than ~ 10^ yr. 

3. On the basis of these results, we predict the observ- 
able properties and the evolution of QSOs at z > 6 
as follows; (i) The evolution of the QSO luminosity 
strongly depends on the redshift Zi at which the ac- 
cretion onto a seed BH got started. In other words, 
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the bright QSOs given z (e.g., z = 6) favor the late 
growth of SMBHs (i.e., z\ « 10) rather than the 
early growth (i.e., Zj « 30) because the timcscalc 
of gas supplied from hosts is shorter than that of 
the star formation in the CND. (ii) In the QSO 
phase, after which the gas supply stopped and when 
the AGN luminosity is above the threshold one, 
^AGN.th (e.g., ^AGN.th = 10^'^ergs~i), there exists 
a stellar rich massive CND with Mg/Mdyn ~ 0.1 
where Mdyn is the dynamical mass of the CND 
plus BH system. On the other hand, in the proto- 
QSO phase (the early phase of a growing BH), 
the proto-QSO has a gas rich massive CNDs with 

Mg/Mdy„ « 1. 

4. Our theoretical model predicts that the distant 
QSO SDSS J1148+5251 discovered at z = 6.42 fa- 
vors the model of late growth of SMBH with Zi = 8 
rather than early growth with Zi = 25. If this is 
the case, the QSO harbors a massive CNDs with 
Mg = 5 X VS^'^Mq and the luminous nuclear star- 
burst in the infrared band with LgB ~ lO^^ergs"^. 
These predictions can be checked by the future ob- 



servations (e.g., ALMA, SPICA, and JWST), in or- 
der to judge whether our model is reasonable for 
the formation of SMBHs. In addition, we predict 
that SDSS J1148-I-5251 have been evolved in a mas- 
sive dark matter halo and/or a massive host (i.e., 
Mdm ~ 1O12-13M0 and/or Mhost ~ IO^Mq), in 
order to build a SMBH with > lO^M© only by 
gas accretion. The progenitor of SDSS J1148-I-5251 
can be super-Eddington objects and observed as 
the high-z ULIRGs such as SMGs with the in- 
frared luminosity Ljr ~ 10*''ergs~^. The pro- 
genitor has a gas rich CND with the mass ratio 
Mg/Mdyn « 1 and Afg > lO^M©, which will be 
observed by ALMA. 
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